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We calculate using perturbative calculations and Ward identities the basic parameters of the 
Fermi Liquid: the scattering vertex, the Landau interaction function, the effective mass, specific 
heat, and physical susceptibilities for a model of two-dimensional (2D) fermions with a short ranged 
interaction at non-zero temperature. The leading temperature dependence of the spin components 
of the scattering vertex, the Landau function, and the spin susceptibility is found to be linear. 
We find that the standard T = relationships for a Galilean-invariant Fermi Liquid are violated 
by finite-temperature terms. The coefficients in the temperature corrections to these relationships 
involve a subtle interplay between contributions from small and large (~ 2kF) momentum processes. 
A connection with previous studies of the 2D Fermi-Liquid parameters is discussed. We conclude 
that the linear leading temperature dependence of the parameters is a generic feature of the 2D 
Fermi Liquid. 
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I. INTRODUCTION 

For more than a decade two-dimensional (2D) fermion systems have been a subject of great interest. One important 
motivation has been the non- Fermi-Liquid behavior observed in high- Tc superconductors above TcB In this context 
the existence ajid stabihty of the Fermi Liquid (FL) in d > 1 has been extensively investigated within various modern 
approaches It seems fair to say that no compelling theoretical evidence has been found for non- Fermi-Liquid 
behavior in systems with short raiiefid interactions and in the absence of coupling to gauge fields. The FL is stable, 
provided that standard conditionsoEl are satisfied. 

Given this a natural Question arises: suppose 2D fermions are in the FL phase in the sense of the Landau's Fermi 
Liquid Theory (FLT) uu We understand "zero order" FL quantities as those defined in the limit of zero temperature 
and zero momentum-energy transfers. Then what is the low-energy behavior of such quantities? In this paper we 
study the leading temperature corrections to the zero-order behavior of the the components of the Landau function 
and the scattering vertex, the quasiparticle's effective mass, and the uniform response functions (compressibility, spin 
susceptibility). Ln-Renormalization Group (RG) language we are studying the corrections coming from the leading 
irrelevant terms. B'tiiy 

Interest in the issue grown in recent years also because these leading corrections providc-thc "bare" temperature 
dependence of the para me | tej rs in theories describing quantum critical phenomena in metalsO In particular, puzzling 
data in several materialal3il3 led to explanations based on an unusual underlying temperature, momentum or frequency 
dependence of electronic susceptibilities.EiO 

The leading temperature corrections to the parameters of a stable FL been, studied for many years. Rather surpris- 
ingly, the issue remains a subject of controversy. For example, it was foundll3 that the leading t epipa rature correction 
to the specific heat coefficient 7 = C/T was T^lnT in d=3 spatial dimensions and T in d=2.tZ[lla. Recently these 
results were also rederived by multidimensional bosonization (see, e.g., Rcf. [ D). Whether the spin and charge sus- 
ceptibilities display similarly anomalous (i.e., non-T^) temperature dependence is the subject of a lengthy controversy 
in the literature: see, e.g., Ref. [ discussion and references therein. FO|i^ most recent reassesment of such results 
see Ref. [ ^ . The prevailing conclusion was that of Carneiro and Pethickij who found no leading In T correction 
to the spin susceptibility of the 3D FL. Their analysis implies that terms cx T are absent in 2D. 

The heuristic argument which runs commonly through the literature to account a posteriori for the absence of 
anomalous terms (in T or in q) in response functions is that although these terms are known to occur in the individual 
self-energy and vertex diagrams, they cancel in final results for symmetry reasons mathematically expressed via Ward 
identities. As we now discuss, this argument is misleading. Although exact symmetries of a model (e.g., global 
gauge and/or rotational invariance) result in relationships between the self-energy and the vertex through the Ward 
identities (cf. Sec. |^ below), these are usually not enough to close the system of equations and demonstrate explicitly 
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the cancellation of the self-energy and vertex corrections. Additional symmetries like, e.g., the chiral symmetry of the 
ID Tomonaga-Luttinger model, are needed in order to do it.Eil In d > 1 analogous extra asymptotic symmetries first 
noted by HaldanecS occur in the fermionic Wilsonian low-energy effective action under some model's restrictions. In 
the recent series of papers by Metzner, Castellani and Di Castro (reviewed in Ref. [ ^) the relationship between the 
extra asymptotic symmetry of the effective action and the RPA (FLT) results for the response functions was clearly 
established. They also clarified a close connection between this approach and multidimensional bosonization.Q The 
additional Ward identities are properties of a model in which all processes involving momentum transfers greater than 
a cutoff 

Aca < kp (1) 

are discarded. These restrictions appear naturally in the Wilsonian low-energy effective action with the UV cutoff 
(|l|) obtained from an RG approach such as that of Shankarp in which interacting fermions are treated by progressive 
elimination of modes towards the Fermi surface. Only asymptotically (i.e., in the limit Aca/kp — > 0, T pjiO) this 
forward-scattering action possesses extra symmetries U{1)°° (and SU{2)°° for rotationally invariant case).E3'cl (It is 
assumed that the initial microscopic action (Hamiltonian) has only ordinary U{1) [and SU{2)] symmetry.) However, 
these extra identities do not constrain possible T- or q- terms in physical quantities, coming from irrelevant terms. 
Even in ID, e.g., the SU{2)l ^ S'C/(2)fl-breaking marginally irrelevant term coming from backscattering results in 
1/ In T- leading correction to the uniform spin susceptibility of the Luttinger liquid.Lj In higher dimensions there are 
always [even in the limit ([^)] special configurations of momenta near the Fermi surface .which give rise to irrelevant 
terms in the effective action which strongly interfere with the forward scattering channelij and invalidate asymptotic 
symmetries. In the framework of the low-energy effective action approach there is also another potential source of 
temperature corrections, i.e., the possibility of T-dependence of the Wilsonian action's vertices (couplings) developing 
on the previous stages of modes elimination before one reaches the effective action scale (|l|) . Apparently, this question 
has not been carefully studied. 

We note that in the most of existing literature studying the leading FL corrections it is assumed that the crucial 
coupling is between quasiparticles and long- wavelength collective modes, i.e., only small momentum scales are taken 
into account. However the possibility of "2fci? singularities", i.e., anomalous temperature terms coming from procesSjes, 
involving large (~ 2kF) momentum transfers, has been pointed out by Misawa for 3D FL already in early TOs.cil 
Apparently due to the lack of experimental evidence of a In T term in the susceptibility of a generic 3D FL and 
also because Misawa's results rely on the analysis of selected diagrams [cf. previous paragraph], they were widely 
disregarded in favor of those of Carneiro and Pethick. In the context of semiconductor physics Stern was the first to 
noteES that in a 2D electron gas the electron scattering rate was proportional to T due to 2kp effects. The consequences 
of the 2kp effects for the leading T-dependence of 2D FL quantities have not been considered in the literature until 
recently. 

TheJssue of the leading correction to FLT has recently been revived by several papers. Belitz, Kirkpatrick and 
VojtaEj presented perturbative calculations, mode-coupling arguments and power counting estimates which showed 
that the leading q dependence of the spin susceptibility (but not the charge susceptibility) was \q\ in 2D (g^lng in 
3D). They did not find the analogous T-correction explicitly, but concluded that one should generally expect a linear 
T-term in the 2D FL susceptibility (T^JirT in 3D). This dependence has important implications for the theory of the 
quantum critical metalliC|Jerromagnet.B3 

Senechal and one of usli3 predicted the occurence of the linear T-corrcctions to the FL vertices from one-loop RG 
calculations based on a 2D effective action. Those RG calculations of the irrelevant corrections rely more on the 
low-energy effective action's phase space constrains in their way to sort out the effective interactions and the scales 
involved, rather than on doing it according to the strength of the latter, as a perturbation theory does. That makes 
desirable to easily obtain a perturbative signal of the RG predictions. Also, an analogous RG calculation of other FL 
quantities was not done. |— . 

Hirashima and TakahashiE3-|performed numerical analyses of perturbative expressions which appeared to confirm 
the prediction of Belitz et oEj for 2D susceptibility. However due to numerical difficulties in handling divergences 
in some terms they were unable even to determine the sign of the coefficient in the leading g-term. Also contrary 
to Belitz et al who focused on the long-wavelength contributions, the authors of Ref. [ ^ emphasized the crucial 
role oij2kp contributions in their findings. Following this Misawa conjectured a phenomenological form for the free 
energyta which results in the linear T-terroJ-n the 2D spin susceptibility and in the coefficient 7, and agrees well with 
the numerical calculations in lowest ordered 

This previous work has left a number of important questions unresolved, including analytic treatment of anomalous 
terms (importarrt for verifying cancellations between different contributions), the relationship to FLT and to Ward 
identities, and the connection to the extensive literature from the 1960-70s. To elucidate the issues in the most trans- 
parent way, we apply the perturbation theory for 2D contact-interacting spin-i fermions, starting from a microscopic 
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action. Although the Landau FLT is not a perturbative theory, for sufficiently weak interactions (assuming the inter- 
action being repulsive and we are above the Kohn-Luttinger temperature) one should be able to find the parameters 
of the FL in terms of the coupling series. 

We present what is apparently the first analytic calculation of the leading T-dcpcndence of the effective mass, 
Landau parameters and response functions of a 2D electron gas, to second order in the interaction strength, including 
all channels and all momentum processes (scales). We take into account the Ward identities explicitly. As we show 
in this paper, the processes involving large 2A;f) momentum transfers are crucial to the anomalous temperature 
dependence of FL quantities. 

The rest of this paper is organized as follows. Section [I A defines the model. Section II B defines the four-point 
vertex to be calculated and its relationship to the Fermi Liquid parameters. We also give there the basic equation 
for that vertex in the one-loop approximation. In Section III we present and discuss our results for the scattering 
amplitude and the Landau function. In Section IV we give the Ward identities which are used for the following 
calculations. In Section^ we calculate the effective mass. Sections VI, VII present the results for t he co mpressibility 
and the spin susceptibility, respectively. The results are recapitulated in the concluding Section VIII . Appendices 
A and B contain detailed presentation of the calculations of the one-loop contributions entering the equation for 
the vertex. One of the goals o f prese nting these technicalities was to give the idea of how more involved two- and 
three-loop calculations of Sees. |V-VII were carried out. I— , 

Some of the results presented here were announced in a previous short communication.Ea 



II. MODEL AND FORMALISM 



A. Model 

We treat interacting fermions at finite temperature in the standard path integral Grassmannian formalism. The 
partition function is 

2 = J Vti^V^P e^o+^'"* (2) 

where the free part of the action is 

^0= / M^)[iuJi+fi-e{ki)]Ml) (3) 
J(i) 



and 

1 r dk, 



(i) PJ i2nr 

(i) = (k„c^,) = /C, (4b) 

where (3 is the inverse temperature, /i the chemical potential, uJi the fermion Matsubara frequencies and ipaii) a 
A'^-component Grassmann field with a spin (flavor, if A'^ ^ 2) index a. Summation over repeated indices is implicit 
throughout this paper. We set fc^ = 1 and h= 1. The S't/(A^)-invariant quartic interaction is 

Sint = -] [ Vic.(l)#(2)V7(3)V'e(4)f/^f (l,2;3,4)5('^+i)(l + 2-3-4) (5) 

4 J(l,2,3,4) 

Here the conservation of energy and momentum is enforced by the symbolic delta function 

+ 2 - 3 - 4) = /3{2T:fSiki + ka - kg - k4)<5^,+^,-^3-^„o (6) 

The SU{N) extension of the physical SU{2) symmetry is useful for different applications,! and the case of the spin-i 
electrons is recovered by setting N — 2. We can decompose the interaction by factorizing its symmetric and antisym- 
metric frequency-momentum- axid spin- dependent parts in the way which explicitly manifests the antisymmetry of 
the function U under exchanged 

t/"/ = f/^/^i' + U^T^f (7) 



3 



where the functions and have the symmetry properties 

[7^(1, 2; 3, 4) = -C/^(2, 1; 3, 4) = -C/^(l, 2; 4, 3) (8a) 
?7^(1,2;3,4) = C/^(2,l;3,4) = C/^(l,2;4,3) (8b) 

while two operators / and T, which are respectively symmetric and antisymmetric in the spin space, are defined as 
follows: 

/"f = 5aeSi3j + 5a^5i3e (9a) 
T^f = S^eSp^ - Sc^Spe (9b) 

If we assume an instantaneous point interaction between electrons, then — and = C/q- The sign of the 
interaction is chosen such that [/o > corresponds to repulsion. Then 

SLt^-\uJ ^„(l)^^(2)V^(3)^,(4)T^f<5('*+i)(l + 2-3-4) (10) 

4 J(l,2,3,4) 

In this paper we consider mainly 2D electrons with the bare spectrum of a free gas e(k) = fc^/2m and the circular 
Fermi surface, interacting via (weak) contact repulsion (^0|), but discuss the consequences of generic spectra and a 
non-circular Fermi surface. 



B. Four-point IPI vertex and parameters of Fermi Liquid 

We define the two-particle Green's function as: 

G2(l,2;3,4) = -(^„(l)V'/3 (2)^^7(3)^^(4)) (11) 
The IPI vertex F is related to G2 as 

G2(l,2;3,4) = (^„(1)^^(3))(V'^(2)4(4)) - (^c.(l)4(4))(^/3(2)^^(3)) 

(V. (l)V5a' (1')) (#(2)^;3' (2'))r^3' (l'> 2'; 3', 4') (^y i^')^,^ i^e' (4')^e(4)) 

l'.2'3'4' 

x5('^+i)(i' + 2'-3'-4') (12) 

where due to the symmetries of the model we can write the one-particle Green's function as 

- (^„(1)V;M2)) = G{\)5^p6^''+'\l - 2) (13) 

Perturbatively, the vertex is constructed by considering only connected one-particle-irreduciblc (IPI) diagrams with 
amputated external legs. In lowest order F^*^^ (1, 2; Z,4) — JJ (1, 2; 3, 4). To shorten notations we use the hat meaning 
that it is an operator in the spin space, and it also comprises two components. Along with the representation (|7|) 
with the components (F'^, F'^) for the vertex, which explicitly shows its antisymmetry, we will also use another 
representation more convenient for some applications, namely separating the vertex into charge and spin components 
via 

f ^ F^^f = r^/,t -f F^r;/ = -^v^^s^.s^, - ^r^pA^^A^, (w) 

where A° (a = 1, iV^ — 1) are Hermitian traceless generators of the SU (N) group, coinciding with the Pauli matrices 
for the case N = 2, and normalized such that 

tr(rA^) = 25"'' (15) 

The components of the vertex in different representations are related as 

Pch = (A^- l)r'^ - (iV+ 1)F'^ (16a) 
Fsp = -F^ - F-^ (16b) 
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Taking into account momentum and energy conservation, we write the vertex as: 

f(l,2;l + Q,2-Q) = f(l,2;Q) , 

with the transfer vector 

Q = 3 - 1 = (q, 



(17) 



(18) 



where «f2 is a bosonic Matsubara frequency. In what follows we will fr'cquently use the property of the components 
Fch and Fgp which is a consequence of the vertex's antisymmetry 



r"f(l,2;Q)=rf"(2,l;-Q) 



reh(l,2;Q) =reh(2,l;-Q), rsp(l, 2; Q) = rsp(2, 1; -Q) 



(19) 



For the calculation of physical quantities of the FL we need f (1, 2; Q) at Q = 0. As is well known0~§l the limit Q ^ 
is not unique, since the vertex is non-analytic function at Q = 0. Two vertices 



F«(l,2) = lim F(l,2; Q) 



(20) 



F"(l,2) = lim 



r(l,2;Q) 



9=0. 



(21) 



can be defined unambiguously at 1 ^ 2 and then continued to 1 ^ 2. (See, e.g., Refs. [ p|,p^ on this subtlety). 
In the calculation of the FLT vertices (real electrons, N = 2) the external momenta are chosen to lie on the Fermi 
surface, and the external Fermionic frequencies are put equal to the minimal Matsubara frequency inT. In the case 
of a rotationally invariant Fermi surface the vertex dependence on the external momenta ki and k2 lying on the 
Fermi surface can be parameterized by the relative angle 9i2 between those momenta and then r'?(6'i2), T^{0i2) can 
be identified with the scattering vertex and the Landau function, respectively. Namely, for the components of the 
scattering vertex {A, B) we have (cf., e.g., Ref.[P]) 



(22) 



SdKp "^m* / {2ttY is the d-dimensional renormalized electron density of states per spin on the Fermi level. 



where Vfj 

Sd is the d-dimensional area of the unit sphere (for d 
So 



2: V]^ = m*/27r), and Z is the field renormalization constant. 



(23) 



B{e,2) - vi,z^Ti^^{e,2 



(24) 



Two components (i^, G) of the Landau interaction function are defined by analogous equations, with the substitution 
q^n,A^F,B^G. 

The one-loop approximation F^^-* of the straightforward perturbation theory for r(l, 2; Q) in diagrammatic form is 
given in Fig. 1. For the model given by (|^,^) we obtain after some spin summations 



fW(l,2;Q) = -^L_C/2/- 



with 



C(l + 2) C/2 \t 



(25) 



L± = L(Q)±L(Q') 



(26) 



We denote Q' = 2 — 1 — Q. The functions L and C coming from the contributions of bubbles (ZS, ZS') and (BCS), 
respectively, are defined and calculated in Appendices ^ and ^ From the one-loop approximation (^) for the IPI 
vertex F valid for generic momenta, energies and transfer, we can obtain now the Fermi Liquid parameters (vertices) 
using the definitions (|2C|,|2l|,p|,p3) . 
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III. TEMPERATURE DEPENDENCE OF THE SCATTERING VERTEX AND THE LANDAU 

FUNCTION 



Taking the appropriate zero transfer limit we obtain from Eqs. ( p5| , p3| , p4| ) the components of the scattering vertex 

A{0)^u~u^ + 2x{q'W-C{qs)u\ (27) 



B{9) = -11-11^ + C{qs)i 



(28) 



with q' ~ 2fci? ainO and qs — 2kp coa9. We absorbed the factor j/q = m/27r into the definition of couphng via 

u = lyoUo . (29) 

Definitions and detailed calculations of the functions involved in the equations above are given in Appendices ^ and 
p. For the components of the Landau function we obtain in the same way 



F{e)^u + 2x{q'y -Ciqs)u^ 



(30) 



G{0)^-u + Ciqs)u^ (31) 

Note that while calculating the FL vertices ( p3|j24| ) at the one- loop level we can put vji = (i.e., m* = m), Z — 1, 
and since L{^1, 0) = 0, the ZS contribution to F, G is zero in the Hmit (|l|). 

Using the results ( A27, A3C ,B3C ,B31) of the Appendices we obtain for the Fourier components of the two vertices: 



Ao = M-u2(ln2A-l + ^r2 



(32) 



Fn = u - it^ In 2A - 2 



(' 



96^ 



(33) 



So = -M + ( In 2A - 1 - 



(34) 



In 2. 



Go = -u + V? In 2A - ^To 



(35) 



V2 96 ° 



(36) 



ln2. 



-To 



(37) 



The terms omitted in the above results are 0{u^Tq) at most. Within one- loop accuracy the standard RPA-type 
relationships between the FL parameters An — Fn/{1 + F„), i3„ — G„/(l + G„) hold. 

Note that taken separately, each contribution of the ZS'- or BCS bubble ha s a leading linear tem perature-dependent 
term in the first two Fourier components of the vertices [cf. Eqs.(A27, A3C) and Eqs.(B3C ,B31), resp.]. However, a 
cancellation of such terms coming from two graphs occurs in the "charge sector" (i.e., in A, F components), while 
the linear T-terms survive in the "spin sector" {B, G components). 

To understand where the temperature dependence of the vertices comes from, note that x(9') ^ ^ function of q' has 
s ome "temperature smearing" near q' 2kp, otherwise it is virtually indistinguishable from its zero-temperature limit 
( [A13| ), provided the condition ( A25 ) applies. So the temperature dependence of the ZS' (exchange) contribution to the 
Fourier components of the vertices comes from integrating of the function x around the "effective transfer" q' ~ 2kF 
in the ZS' graph, i.e., when incoming momenta ki ~ — k2, while |ki| — |k2| — kp. Notice that the transfer q in the 
proper sense (^Jl8) is zero. This source of temper ature dependence was not considered by Pethick and Carneiro, 
for example. In the same vein, one can see from Eq.( |B15 ) that the temperature dependence of the BCS contribution 
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comes from regions of small Qs, i.e., again when ki ^ ~J^ 2- T he co efBcients in the temperature expansions of the 

Fourier components of x and C are such [cf. Eqs.( A27 A30 , B30 B31 )] that the cancellation occurs only for the linear 
terms in the charge components of the vertices, while the higher order terms in the temperature survive. Since to the 
best of our knowledge there is no particular reason for this cancellation of the linear terms to happen, we expect that 
linear T-term of the charge component (s) would come out of a more realistic calculations of the FL parameters. 

Let us elaborate on the last point a bit more. First notice that in the perturbative approximation (^) for the 
vertex in the model with a constant couplings Uq (bare vertex), all the three one-loop terms have the same factor 
Uq in front of the bubble contributions L, C. Had these coupling factors been different for each of the three graph, 
the cancellation of the temperature corrections would not have happened. On the other hand, even the simplest 
approximation ( p5|) shows that fermionic interactions drive any bare coupling constant towards a coupling function 
of, generally speaking, three momenta [cf. Eq.(p7|)]. More advanced calculations of the FL parameters, like, e.g., RG, 
would reproduce such development of the momentum dependence in couplings in the regime of effective action (|^). 
This would prevent the cancellations between different terms, e.g™ between the one-loop ZS' and BCS contributions, 
which we have obtained in the simple perturbational framework.Eil 

Another artifact of the naive perturbative calculation of the vertex is its ultraviolet divergence coming from the 
BCS bubble. The In A-term makes our perturbation theory applicable only if the interaction is small enough that 
the arbitrary, generally speaking, ultaviolet cutoff A can be chosen to satisfy 1 <C A < ^e^^", (u < 1). However, we 
should not worry too much about this issue, since more sophisticated calculations based, e.g., on summations of the 
diagrammatic series or an RG would result in an effective screening of the repulsive interaction in the BCS channel. 
See, e.g., Refs.[ ^,^,||. (Although providing the possibility of the infrared instability (divergence) ^exponentially low 
temperatures (much lower than, say, we are working at) through the Kohn-Luttinger mechanism.Ej) 

The linear dependence of the leading temperature corrections to the FL parameters seems not to be sensitive to the 
actual form of the one-particle spectrum. The same dependence (apart from model-sensitive prefactors) was obtaiiied 
in the previous RG analysis of the effective action model for 2D spinless fermions with a linearized spectrum.E2l In 
that study the integration of the "effective" transfer in the exchange (ZS') graph contribution to the RG flow was 
also resulting in the linear T-term in the vertex corrections. According to Ref. [ such temperature behavior can 
be understood from dimensional arguments. 



IV. WARD IDENTITIES 



In order to derive the Ward identities, we icdlowed closely the methodology of Ref.[ adapted to the path-integral 
technique according to standard approaches. The partition function is given by the path integral (|^). In coordinate- 
imaginary-time space the free part of the action is 



5*0 - / dXij^iX) 



d A 

OT 2m 



MX) (38) 



We consider pair interactions which preserve the total spin and number of particles 

Sint= f MXl)MXl)U{^i-^2)S{Ti^T2)MX2)MX2) (39) 

■J{XuX2) 

Here J^^-^ = dr J dx. Notice that interaction (^o|) is a special case of (|39|). Then the following charge Ward identity 

.l^-^(2k,-q) ='-N J,,, .17-^(2k,-q)^"/^(^' ^' + 2) (40) 

and the spin Ward identity 



can be derived. From (|14[|15|) one can find 



K^f) = -^reh (42a) 

A»,r^;A^„ = -2{N^ - i)r,p (42b) 
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These equations are helpful in order to unclutter the r.h.s of the Ward identities. A useful consequence of the Ward 
identities is the the identity: 



[reh(l, /C; Q) - r,p(l,/C; Q)]G(/C)G(/C - Q)- 



r(2k-q) 



zr!-^(2ki+q) 



= 



[r^(l, /C; Q) - r^(l, /C; Q)] G(1C)G{1C - Q)- ^q"',;. = 

{K) «f^-^(2ki+q) 



(43a) 
(43b) 



The above relationshiKi are non-trivial, since rch(l,/C; Q) ^ r,p(l,/C; Q) [r^(l,/C; Q) ^ r'4(l,/C; Q)]. There is also 
another Ward identityO 



which follows from the Galilean invariance. 



k2rch(l, 2; -in)G{2)G(2 + iVl) 



(2) 



V. EFFECTIVE MASS 



A. FLT calculation 



(44) 



In this subsection we calculate the effective mass using the Ward identities. The self-energy is defined such that 

G-i(l) = GoHl)-S(l) (45) 

From now on we assume Go(l) in the above equation to be not the exact bare Green's function of the non- interacting 
system, but that with the shifted chemical potential. We gauge the shift to be such that it removes all tadpole 
insertions to the Green's functions in diagrammatics. 
The effective mass to* is given by the following equation 



TO 

m 



1 



din 



1 



c?k kp 



(46) 



Since we are working with the finite-temperature Matsubara Green's functions and loi is a Fermionic Matsubara 
frequency, we understand ioji in the sense that the zero-frequency limit will be taken after the appropriate 
analytical continuation of the final result for to*(1). In the second order of the perturbation theory only the "sunrise" 
self-energy diagram contributes to the mass renormalization. In this approximation 



TO 
TO 



1 



(9I](1) TO aE(l) ki 
diO. kp 9k kp 



ki sSf 



The charge Ward identity can be written as 

'(2) 



ak 



/ i^reh(l,2;-q)G(2)G(2 + q) 
J{2) m 



while the Galilean-invariance Ward identity (|4^) reads 

as(i) f kik2 



din 



rch(l,2;-if2)G(2)G(2 + ifl) 



'(2) ft-l 

Combining the three above equations together, we get 
kik2 



TO 
TO 



1 



(2) f^F 



rch(l,2;-il7)G(2)G(2 + irj) 



rch(l,2;-q)G(2)G(2 + q) 



(47) 



(48) 



(49) 



(50) 



Since the above equation has the accuracy 0{U^), we can use the one-loop approximation for the vertex, while the 
Green's functions can be approximated by Gq. 
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For the the contact interaction (|10|) {N = 2) the vertex's one-loop approximation ( |25| ) gives 

reh(l, 2;~Q) = Uo + L{Q)Ul - 2L(1 - 2 - Q)U^ ~ C(l + 2)Ul (51) 

It is easy to verify that only the third and the fourth terms in the above expression (coming from ZS' and BCS graphs, 
resp.) give non-zero contributions to m* . Plugging in the formulas for the ZS' and BCS graph contributions and 
carrying out some manipulations, we obtain (recall that iloi ^0, ki e Sp) 



m 



1 + -u^ 

m 



1^^2 f kik2 rtanh(i/3^k3) - tanh(i/3,^k3+k2-ki) tanh(i/3^k3) 1 P 



j- tannt^^P^ka; - tanni_^p,^k3+k2-ki j tanni^^/^ska j 1 

4 ° i(k2,k3) -iWl + ^k2 + ?k3 - '?k3-|-k2-kl -JWl -^k2 +Ck3 +Ck3-k2-kl J COSh^(i/3^k2) 

(52) 

Note that without applying the Ward identities, the above equation can be derived directly from Eq. (^) where the 
self-energy is given by the sunrise diagram. However, it takes much more tedious and lengthy calculations to verify 
the cancellation of other terms appearing at intermediate steps of the analysis. 

Before proceeding further with the calculations, let us first comment on result (|5^). In the the zero-temperature 

limit /?/ cosh^(|/?^k2) oc '5(^k2) enforces k2 to lie on the Fermi surface, while the integration over the angle kik2 is 
equivalent to the calculation of the first Fourier harmonic of the expression in the curly brackets. The latter we have 
calculated already, and it is the sum of the ZS' and BCS contributions to the FL vertex. (Recall that the first order 
term and the ZS graph do not contribute to the first Fourier componeiii of the vertex.) Then Eq. (|5^) in the the 
zero-temperature limit readily recovers the standard result of the FLTUIj 

"^*^^ = °^ =l + fi(T^O) (53) 
m 

A straightforward extension of Eq. ( ^3|) to finite temperatures like m*(T)/m = 1 + Fi{T) is not valid, since according 
to Eq. ( p2| ) ■m*{T) contains an extra contribution from the "off-shell" integration over ^2 (Ck2) normal to the Fermi 
surface, albeit the factor /?/ cosh^(i/3^k2) niakes this contribution well-localized near the Fermi surface. In other 
words the vertex entering the r.h.s. of Eq. (p2h is not exactly the FL vertex F(T) as it is defined in the FLT (cf. 



definitions in Sec. HE), since one of its momenta (namely, k2) is not confined to the Fermij-surface. 

The calculation of the temperature expansion of the effective mass m* is quite involvedHJ We did not go beyond 
the linear-temperature terms. Evaluation of the contributions of two graphs to m* denoted as 



gives 



— = l + u'{Mzs' +Mbcs) (54) 
m 



Mzs' ^2il + \n2)To + 0{T§) (55a) 
Mbcs = ^ - J(l + ln2)To + 0(e-i/^«) (55b) 



Thus we get 

m* 1 



^ = 1 + 2^ + ^("^0 ) (56) 



(Recall that dimensionless coupling and temperature are introduced by (29) and ( A25 ), resp.) In a close analogy 
with the cancellation of the linear-temperature terms in the Fourier components of the FL vertices A and here the 
cancellation occurs between additive linear-T corrections coming from both "on-shell" (i.e., linear T-term coming from 
the 2fcF-contribution to the vertex) and "off-shell" (i.e., small- momentum contribution) integrations in two diagrams. 
Moreover, the "on-shell" ("oflt-shell") T-term of the ZS' graph cancels the "on-shell" ("off-shell") T-term of the BCS 
graph, correspondingly. 



B. Alternative Calculation 



In this subection we present an alternative direct evaluation of the leading temperature correction. The advantage 
of this evaluation is that it applies also to non-Galilean invariant situations and helps to more clearly establish the 
relation to previous studies. It also indicates possible generafizations. 
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We begin from the explicit expression for the self-energy correction given by the "sunrise" diagram 

I](P) = [/2/" G{V + Q + ]C)G{-]C)G{Q) (57) 

This choice of variables is convenient because the T— linear contributions to the effective mass will be seen to arise 
either from regions where q+k is small or from regions where it has magnitude 2pp. The bare electron Green function 
G is given by ([A^). Evaluation of the Matusbara sums leads to 

^^p^,^)^U'[ K(6)-no(^,)][6fe-6)+»o(W0] (58) 

Jk,q + iq~ S.k- Cp+k+q 

and by using [no{S,k) — no{^q)]\b{S_q — ^j.)] = no(fq)[l — ?io(Cfc)] ("-o and b are the Fermi-Dirac and Bose-Einstein 
distribution functions, rcsp.), and relabelling variables we obtain 

np,^^) = Ul^ ( no(Oo(a)j) (59) 

with J = 2/{iuj + £.q-£,k- Cp+k+q) + l/{-iuJ + (.q+ 6 - Cp+k+q). 

Now the first term in the expression for S cannot give rise to a term of order T because the integral of is not 
confined to the region near the Fermi surface. Indeed, differentiating this term on temperature leads to a numerator 
containing the factor ^q. Because the integral is not confined to the region near the Fermi surface, obtaining a 
nonzero value for the ^q integral would require an additional factor of ^q divided by a scale of the order of Ep, meaning 
that the integral would be of order as in the usual Sommerfeld expansion. 

We therefore focus on possible T-linear contributions from the second term. To obtain these we will make the 
assumption, to be justified a posteriori, that all energies {^q,£,k,£,p) are near zero. We parameterize the q, k integrals 
by £,k,q and angles. We choose 9i as the angle between p and q + k and 02 as the angle between q and k. For a 
spherical Fermi surface we have 

ep+k+q = A + Bcos(0i) (60) 

with, (up to first order in deviations from the Fermi surface) 

A =[AEF+^q+^k] cos'' (02/2) (61) 
B = [AEF + 2^p + ^q + ^k]cosie2/2) (62) 

We may now do the integral over 9i obtaining 

f 

S(p,«w) = -— d^qd^kd92no{^k)no{^q)K (63) 

with 

^ ^ 2sgn{Ci) ^ sgn{C2) 

and 

Ci = ^ + - a - Cp - [^Ef +£.q+ a] cos2((?2/2) (65) 
C2 = -^ + ^q+ik-£.p- [^Ef + + a] cos2(02/2) (66) 

The crucial point is now that (after analytic continuation in lo) the region of phase space in which the square root is 
real is very small, and in particular corresponds to 6*2 ~ or 62 ~ ±7r, as well as to small values of the energies, thus 
justifying the assumption made above. Let us consider these two regions separately. 

• 62 ~ Ti'. These angles correspond to processes,ia which a small total momentum is transferred to the electron, 
i.e. to the processes which are usually consideredlHHlSI to give rise to anomalous terms in C/T, etc. In this regime we 
write O2 = ±7r — y; to acount for ±7r, y may have either sign. To leading order we have 

Ci - CO + - a - (67) 
C2 = -tJ + Cg + a - ip (68) 
B = 2EFy (69) 
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Note that if all the ^ are small, then 6*2 is confined to small angles, where the approximation is accurate. Performing 
the angle integral gives 



'(P,^) = / d^,d^kn.{CMQL^"^-^^ (70) 



with 



j^small ^ 2.sgn{u + " 6 " Cp) + Sgni-CO + + 6 " Cp) (71) 

The remainder of the evaluation is simple: to obtain the quasiparticle velocity renormalization we must take the sum 
of the w and derivatives; thus the first term does not contribute at all while the second term gives — 2(5(^q + S.k)- 
for the Lu derivative and the same for the S,p derivative. Substituting this into the integrals and evaluating gives 

\srnall = ^ / d&"-o (Cfe )"o (-& ) = (72) 

• 02 ~ 0: These angles correspond to processes in which a momentum near 2pp is tranferred to the electrons. In 
this regime we have sgn{Ci) — sgn{C2) = —1. Further, as both Ci,2 and B are of order Ep we write — = 
{C + B)*{C-B) = -8Ef*{C + B). We have 

Ci+B = uj-^k+^,+Cp + IefOI (73) 
C2+B^-Lo + ^p + ^, + ^k + ^EfOI (74) 
Again the 02 integral may be done, yielding 

r 

I]2f-(p,c.) = -— J d^,d^kn.ia)n,i^,)L^P- (75) 

with 

L^P- = 2e{Lu - a ^p) + 0{-uj + ^p+^g+ a) (76) 

As above, one of the two terms give no contribution; the other has a factor of two but when differentiated gives only 
S not 26. The sign is opposite. Thus this term cancels the previous term, leading to no T-linear term in the effective 
mass enhancement. 

An almost identical calculation of the order term in the free energy shows that there is no T-linear term in the 
specific heat coefficient. 

As we see from these calculations, the coefficients of the linear T-contributions to the effective mass coming from 
small and nearly 2pp processes are model-dependent. Any model modifications, e.g., momentum dependent interaction 
or non-circular Fermi surface, changes these coefficients, eliminating the cancellation and resulting in apperance of 
the linear T-term in the effective mass. The samCpConclusion applies to the specific heat coefficient. This agrees with 
the recent experimental data on liquid "^He films.Ea 

VI. DENSITY-DENSITY CORRELATION FUNCTION 

In this section we consider the density-density correlation function >c(Q), defined as follows: 

6^'+'\0)>c{Q) ^ {p{Q)pHQ)) (77) 

where the density operator p{Q) is 

P(Q)= / V^a(l)^a(l+Q) (78) 
J{1) 



and p(Q) — p{Q) — (p(Q)) stands for the density fluctuation. Using the definition for the IPI four-point vertex ([T^), 
one obtains 
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^(Q) = -n\ f G(1)G(1 -Q)+ f G(1)G(1 - Q)reh(l, 2; -Q)G(2)G(2 + Q) 

J(l,2) 



The correlation function has the properties 

><(q = 0, 0) = 
>^(Q) = >^(-Q) 

To derive (po|a) it is enough to re-write the charge Ward identity as 



(79) 



(80a) 
(80b) 



G(1-S)G(1) 



G(l - Q) - G(l) 



in 



q(2ki-q) 
2m 

q(2k2+q) 



2m 



(2) ifi - 



q(2ki-q) 



^G(1)G(1 - Q)reh(l, 2; -Q)G(2)G(2 + Q) 



(81) 



Using the above expression for the first term in the brackets on the r.h.s. of Eq. (|79|), one can easily establish Eq. (|80p). 
As to the second Eq. (|80|b), it can be proved with the use of Eqs.(p^). 

In order to find the compressibility of the system we need to calculate the density function in the non-trivial (static) 
zero transfer limit, which we denote = >i'(q ^ 0, = 0). To do it in the second order over interaction by a naive 
perturbation theory, one needs to take into account the sunrise diagram for the self-energy in the Dyson equation ( |45| ) 
and to use such approximation for the total Green's function entering the first term in the brackets on the r.h.s. of 
Eq. ([79]). For the second term in that equation one can use the "bare" Green's functions [cf. comment after Eq. (|4^)] 
and the one-loop approximation for the vertex. In diagrammatic language this amounts to the calculation of a set of 
three-loop diagrams, which is very difficult in practice. The first problem is a proliferation of terms after doing each 
intermidiate step of intagsation (summation) for a given three-loop diagram. The second is that some of these terms 
exhibit IR divergencies .t3 It takes additional efforts to assure the final cancellation of those divergencies between 
different terms. As in the effective mass calculation, the use of the Ward identities drastically simplifies the problem. 
[Cf. comment after Eq. (|2|)]. 

In applying the Ward identities for our calculations we will use the approaches developed in the classic derivations 
of FLTi We denote 



G^(l)^G(l)G(l + q) 
G^(1) = G(1)G(1 + *1^) 
Let us also introduce two quantities An and K as 

G^(1)=G^(1)-Ah(1) 



(82a) 
(82b) 

(83) 



Then from Eqs.d7|,||) we have 



where 



Eqs. (MM,) give 



r?ji,2) = r^ji,2)-if(i,2) 



/ G^(l)+ / G^(l)r^,(l,2)G^(2)->^2 

J{1) J{1,2) 



X2 



(1:2) 



Glil)K{l,2)Gli2) 



GlW+l G2(l)r^Ul,2)G^(2)=0 
Subtracting this from Eq. (|85|) and using (^Sh, we obtain 



(84) 



(85) 



(86) 



(87) 
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N 



(2) 



rg,(l,2)G^(2) 



(1,2) 



Afi(l)rS,(l,2)Afl(2)->.2 



Note that from the charge Ward identity 



mi) 

din 



f r2,(l,2)G2(2) 

J(2) 



(89) 



So far we did not use any approximations, and Eq. ( p8| ) is exact. Now let us simpUfy it at the level 0{Uq) for the 
model with interaction (jl^) {N — 2). From the one- loop approximation for the vertex (^) we see that 



(90) 



Also the Green's functions in Eq. ( |86| ) can be approximated by their bare forms, resulting in 

^2^^u^ + Oiu^) (91) 

ZTT 

After the standard analytical continuation to real frequencies, one can find for the parameter Afl(l), 1 = (ki,ti;i) 



Ai?(l) 



1 



5S(1) m aE(l) ki 



du! 



where 



A(i) = -e:i<5(c.i-ek,; 



kp 9k kp 



(3 

^1 = 7; 



A(l) 



2cosh2(i/3CkJ 



(92) 



(93) 



Eq. ( p2| ) is the extension to the finite-temperature of the well-known result Aij(l,T = 0) = Z'^m* /kp6{uJi)S{ki — 
kp) (cf_, e.g^ Ref. [ ||), combined with expansion of Z^m* to order . Up to the terms 0{u ) we have from 



Eqs. 



mm 92) 



k1 

T 



2tt 2tt 



A(l) 



(1) 



aE(l) m 9S(1) ki 



dijj 



kp 9k /ci? 



(1,2) 



A(l)r^Jl,2)A(2) 



(94) 



As we have shown before for the calculation of the effective mass, the first integral on the r.h.s. of Eq. ( p4[ ) can be 
written via the Ward identities as the second term on the r.h.s. of Eq. ( pO| ) (With the only difference that in Eq. ( p^ ) 
(ki,a;i) is arbitrary.) Since two limits of the vertex differ only by a constant term ( pO| ) which disappears after the 

integration over the angle kik2, we can finally write 



(95) 



_2n 

Ji = — 
m 



kik2 



A(l)r2,(l,2)A(2) 



/o 



2tt 



(1,2) ^^1 

/ A(l)rS,(l,2)A(2) 

"'(1,2) 



Eqs. ( |95| , plD for the response function are manifestly free of spurious IR divergences and readily recover the standard 
FLT results at zero temperature. We recall [cf. Eq. (^3])] that the components of the scattering vertex and the Landau 
function are just the appropriate limits of the vertex Fch times the factor vp^Z^ {vpi = m*/27r). The latter reduces 
simply to ni/2ii in our approximation. Then at T = [cf. Sec. ^ for the explanations of T ^ limit, and note also 
that from Eq. (||) = ^_ 0(^3)] ^^^^ ^^^^ (||) 



^?(T = 0) = -(1 + F2 + Fi-Fo) 

TT 

which is nothing but the expansion to 0{it') of 2D FLT result (cf., e.g., Ref. [ 



^FLT(j. ^ 0) ^ 



ml + Fi 
7l + Fo 



(96) 



(97) 
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Using Eq. ( |5T| ) for the vertex entering the first integral term in Eq. (^Sj), we obtain: 

27r^^2 /" kik2 £1(2:2 rtanh(i/3Ck3) - tanh(i/3^k3+k2-ki) tanh(i/3^k3) 



"1 J(ki,k2,k3) kf 4: I -Cki + Ck2 + Cka - ^k3+k2-ki -^ki - ^k2 + Cka + ^k3-k2-ki J 

_ fZS' , fBCS 
= Jl + /l 

where we indicated expHcitly two contributions to the response function coming from the ZS' and BCS one-loop 
corrections of the vertex. For the second integral on the r.h.s. of Eq. ( |95| ) we have 

fo = ^Uo + fo'''+fi^' (99) 

where the first term is the bare vertex contribution to >f', while the two other one- loop vertex corrections are given 
by the same formulas as in Eq. (|9|) with an obvious replacement kik2/fci 1-^ 1. P-, 

We were able to analytically calculate the above terms in the leading order of their temperature dependence. LJ The 
results are 

(100a) 

Ito) (100b) 



fZS' 

Jo 




fBCS 

Jo 


= w2(-ln2A + 


fZS' 
Jl 




fBCS 
Jl 





(100c) 
(lOOd) 

up to the order 0{u^Tq). Note that the BCS contribution to the term /o is ultraviolet divergent and we regularized 
it in the same way as we did it for that contribution to the vertex, i.e., by introducing an UV cutoff A ^ 1. [Cf. 



Sec. Ill .] Applying Eq. (^9|) for the term fi while passing from Eq. (|9J) to Eq. (|95|) may raise some questions from 
a cautious reader, since the ultraviolet regularization we used to calculate the BCS contribution would break the 
Galilean invariance, which is indispensable for the derivation of (^). However /i is UV-convergent, thus can be 
calculated in the limit of an infinite cutoff when the Galilean invariance is preserved. 
Combining the above results together, we obtain 

^9 = _ u _ 1^2 ^ «2 jj^2A) + Oiu^T^) (101) 

TT 2 

Once again, the leading linear T-corrections which can be traced back to the ZS'- and BCS-loop contributions to the 
vertex cancel in each of the terms /o and /i separately. It is useful to keep in mind that albeit we expressed /i and fo 
in terms of the vertex only, these contributions entangle both "purely vertex" corrections and self-energy corrections 
to the response function. The latter are just expressed in terms of the vertex via the Ward identities. 
Let us finally remindo that the compressibility K = ji'^/n^, where n is electron density. 

VII. SPIN-SPIN CORRELATION FUNCTION 

To find the spin susceptibility in the SU{N) formalism, we consider the spin (flavor, if 7^ 2) correlation function 
X°''(Q), defined as follows: 



^(rf+i)(0);^«&(Q) = (5-(Q)5''t(g)) (102) 



where the flavor density operators 



S^iQ)^^ / AS^^„(l)V/3(l + Q) (103a) 
^ ^(1) 

s^HQ) I Kp^o.{i)M^ - Q) (103b) 
^ ^(1) 



and Q is the gyromagnetic ratio. Using Eqs. (|T2,E3) one can find 
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X^'CQ) = -^S^' \ I G(1)G(1 -Q)+ I G(1)G(1 - Q)r,p(l, 2; -Q)G(2)G(2 + Q)] (104) 

Note that in a paramagnetic state the response is the same along all the (iV^ — 1) directions a. Then we follow 
the same steps as in the previous section, with the only difference that instead of the charge Ward identity we use 
the spin identity (pl|). The latter can be casted in the form of Eq. (81) where Fch ^ Tsp. Function x°'°'{.Q) has 
the same properties ( |80| ) as those of the density function. The uniform susceptibility is given by the static limit 
= -)(°-°-{c^ 0, = 0) (no sum over a). The equation we obtain for 2x'/0^ is the same as the r.h.s of Eq. (|8^), but 



- ch sp 



rE,. Then wc use the consequence of the both spin and charge Ward identities [cf. Eq. ([43|)]: 



/ T%il,2)GU2)= f ri(l,2)G2(2) (105) 

J{2) J(2) 

and afterwards proceed as in the previous section. From the one-loop approximation ( p5| ) 

r,p(l, 2; -Q) - -Uo + L{Q)Ul + C(l + 2)Ul (106) 
Thus up to the terms 0{u^) we obtain 

x' = ^{l + u' + h-9o} (107) 



50^^ / A(l)Ffp(l,2)A(2) 

"'(1,2) 



(108) 



Recalling once again what we said after Eq. ( p5[ ) and noting that from Eq. (^5|) Gq = + 0{u^), we easily see that 
in the zero-temperature limit the above equation gives 

^9(T = 0) = ^(1 + G2 + Fi-Go) (109) 
reproducing thus the second-order expansion of the 2D FLT resulti 

Luckily we do not have to do new calculations since 

50 = - /o''''^ (111) 
So we find the spin susceptibility to have a linear leading temperature correction 

X« = ^(1 + u + ^7.2 _ ^2 i^2A + ^^2) + Oiu'T^) (112) 
47r 2 2 



VIII. SUMMARY AND DISCUSSION 



In this paper we have systematically examined the leading temperature corrections to FLT in two spatial dimensions. 
We find for the model of a 2D electron gas with a contact interaction that to order v? the leading T-dcpendcncc of the 
FL parameters in the spin sector is linear in temperature, while for the parameters in the charge sector and for the 
effective mass a cancellation of the leading T-corrections occurs, and their expansions start from the terms quadratic 
in temperature. 

The particularly interesting result we found is the leading linear temperature dependence of the spin susceptibility 

^^l^u^l- (113) 

According to the perturbative calculations of Belitz et the 2D FL susceptibility has a leading linear correction 
|q| a t T — with a positive coefficient which is of the second opder in interaction, i.e., theirp^sult has a structure 



of Eq. (112). This also agrees with the phenomenology of Misawalla and the numerical results E3 
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Our results reveal the crucial importance of 2kp processes contributions into the low-energy parameters. We remind 
that recovering of the angular dependence of the FL vertices in the whole region where the angle varies, involves taking 
into account transfers varying from zero to 2kp. We have seen it from the direct bubbles evaluation, but the argument 
is non-perturbative. To appreciate this let us write the antisymmetry condition for the vertex as 

r^*^ki,k2;q) -±r^*\ki,k2;k2-ki -q) 

So, small- and large-transfer scales are intrinsically related due to the Pauli principle, and any calculation of the low- 
energy corrections should take this into account. Eventually the integration over whole range of the transfer values 
through a given loop results in the temperature corrections not only to the vertex, but to the response functions as 
well, when that loop is a part of a more complex diagram. The big transfer (~ 2fci?) contributions essentially gave 



rise to the spin susceptibility result (113) in the three-loop approximation. 

We find that the relationships known from the classical FLT derivations at T = for the parameters of Galilean- 
invariant FL (e.g, the effective mass, response functions vz components of the Landau function) are violated by finite- 
temperature terms. The coefficients in the temperature corrections to these relationships subtly involve contrubutions 
from small and large (~ 2kp) momentum processes. 

Concerning the argument for the cancellation of anomalous terms in the response functions due to Ward identities: 
We have calculated the vertices at the one- loop level 0(?i^). Through the Ward identities the self-energy corrections 
were taken into account with the same accuracy. There are no more terms of the order 0{u^) to cancel the temperature 



dependence (113). Thus, the linear temperature dependence of susceptibility (or weaker temperature dependence of 
the compressibility) does not contradict the exact Ward identities known to us, moreover in our results for the response 
functions both vertex and self-energy corrections are included on the same footing by using the Ward identities. 

For more realistic models of electrons in (quasi)-2D crystals, i.e., for various tight-binding spectra and fillings, the 
free-gas-like square-root-2fci? singularities (with kp depending on a chosen direction in q-space) are known to exist in 
the Lindhard functions.123 We think this is enough to result in linear T-terms in physical quantities analogous to what 
we found in this study. We argue that the cancellation of the T-terms in some FL parameters is special to second 
order perturbation theory and the model considered, while the leading linear temperature corrections are a generic 
feature of the 2D FL. 

We hope our results may be experimentally tested in real 2D FL sjfstems. For example, a very naive fit of the 
temperature dependence of the spin susceptibility in Sr2Ru04 systemED when it is in the 2D metallic regime (above 
3D crossover temperature) shows that the data are compatible with the form (pSh. We expect our results stimulate 



a more detailed examination of the leading temperature dependences of response functions in 2D systems. 
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APPENDIX A: ZERO SOUND (ZS) AND PEIERLS (ZS') GRAPHS 
1. Finite-temperature Lindhard function 

We define the Lindhard function as 

L(Q) = L(Q,f^„) = i^ / Go(K,c^,„) Go(K + Q,c^,„ + aO (Al) 

It corresponds to the contribution of the particle-hole loop in graph ZS or ZS', depending on the actual value of Q we 
put. For the purpose of notational economy in the Appendices we will reserve capital letters for momentum vectors, 
while small letters correspond to their dimensionless counterparts defined below. The free electron Green's function 
entering the loop is given by 

Go\K,LUn)=iu;n + ^^-e{K) (A2) 

with its spectrum 
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2m 2m 

A straightforward summation over Matsubara frequencies results in 



e(K) = ^, ^, = ^, CK^e(K)-/. (A3) 



^ / ^ "°(^K) -no(eK.c) 
y (27r)^ illn + Ck - Ck+q 

We denote by Uq the Fermi-Dirac distribution function 

We choose 2D polar axis along the transfer vector such that KQ = KQ cosd. Then the Lindhard function can be 
transformed to the expression 



/•oo /'27r 

L{Q,0) = -^J KdKnoi^K) J < 



The angular integration above is understood in the sense of the principal value. We define the dimensionless Lindhard 
function x zero frequency transfer as 

L(Q,0) = -j/ox(Q), ^0 = ^ (A7) 
where i^o is the free 2D electron density of states per spin. Let us also introduce the new (dimensionless) variables 

/3o = ^,.= ^,fc^^ (A8) 
4m ' ^ 2Kf Kf 

In order to proceed we use the following result for the principal value of the integral /q at a > 0: 

e(a2 > 1) (A9) 



/o a? - cos2 i9 aVa^ - 1 
where is the event (Heaviside) function. After the angular integration we obtain 



1 /■■? hrik 

X{q) = / ^=^{1 - tanh/3o(ft2 - 1)} (AlO) 



The above formula can be easily casted into the final form 



,(,) = - + -tanh/?„--|^ -^^1--^ (All) 
To make zero-temperature check of our result the following formula is useful 



/3o-oo: I -^5{z) (A12) 

^ cosh Poz 



Then from Eq. (All) we easily recover the result known for T — 0, which was first found by Stern:E3 

q^<l: x{q)^l (Al3a) 
g2>l: ^{q) = i- ^l-}- (A13b) 
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2. Fourier components 



Let us consider now the Lindhard function at the "effective" transfer 

Q' = K^,') - K^,^) - Q 



(A14) 



which occurs in the ZS' loop when the "real" transfer in the vertex is Q. Parameterizing vectors K^'' and K^'' by 
their angles {0i, 62) relatively to the direction of vector Q, we get 



Q'^ = + 4if| sin^ — + AQKf sin sm 



(A15) 



In terms of new angles 



72 — Cl , 02+ Cl 
7. > ^ = 7. 



(A16) 



and of the dimensionless variables (A8) one reads 

q'^ = + svo? 6 + 2q sin 6 sin 7/) 
Let us consider now xil') zero vertex transfer q — 0, i.e., q'^ = = sin^ 9. For the Fourier transform 

xu = - f decos{2vie)x{e) 



(A17) 



(A18) 



of x(0) given by Eg. ( All ), we have 



X« = ( ^ + ^ tanh/3o )(5„,o - 6xn 



(A19) 



where 



Then 



5Xr. 



27T Jo 



de cos(2n6') 



dz 



_ l + z 

_i cosh^ PqzS siv? 6 



(A20) 



Jo a;v 1 — J-i cosh (ioz j-i cosh fioz Jy 



dx 



Jo x\J\ — 7-1 cosh^ /3( 
Using the result 



cosh f5oz Jym x\/\ ~ X- 



Va;2-(l + z) (A21) 



we get 



— \/ T-^ — = — (1 — a) 



/3o Z"" rf^ 



<5xo = - tanh /3o „ / 2 

2 2 cosh /3oz 



(A22) 



(A23) 



Thus 



Xo 



1 /3o /"^ rf^ 



2 2 7q cosh^ /3oz 
The dimensionless temperature Tq is 



1 I [To dz 



2 2 Jo cosh 



2 Vi - T02: 



(A24) 



1 4mr 



(A25) 
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Taylor-expanding the square root under the integral, we can then extend the upper limit to the infinity with exponential 
accuracy. Thus, 



After simple integration we arrive to the sought result: 

Xo = 1 - i^To - ^T„^ - + 0(r„4) (A27) 

For the next Fourier component using cos 26 — 1 — 2 sin^ 9 we obtain 

J. 2/3o f° dz xdx /-^ ^ (3o zdz . , 

Xi^-Sxo + ^ —J—/ -j==V^^^W+^) = -Sxo + ^ ,2^ (A28) 

TT j_i cosh poz J^/T+^ V 1 " 2: ^ Jo cosh p^z 



Thus, 



2 

and combining the last formula with the previous results, we have 

In 2^ TT^ . 9C(3) , 



XI = -<5xo + + ©(e-i/^") (A29) 



XI = -To - ^To^ - + O(To^) (A30) 



APPENDIX B: BCS GRAPH 



1. BCS loop contribution at finite temperature 

We define the contribution of the particle-particle loop in the BCS graph as 

C(S.) = C(Q„1],) = ^V / -^Go(K,u;,„) Go(-K + Q„-< 

where 

Qs = Ki+K2, ^s^iOi+UJ2 

After summation of the Matzubara frequencies we get 

dK no(CK) -rio(-e„K+Qj _ 1 [ dK tanh(i/5eK)^ 



{2^y 



tanh(i/3^_K+Q, 



(27r) 



$K + ? K+Q, 



(Bl) 



(B2) 



(B3) 



For the case we are interested in (i.e., zero sum frequency fis and vectors Ki, K2 lying on the Fermi surface, cf. also 
notations (A14, AI6)) we have 



K^^^ I = 2A> cos 61, = 



(B4) 



To regularize the momentum integral in (B3) we will introduce an ultraviolet cutoff A which we assume to be the 
largest momentum scale of the prob lem. Under this condition wc shift the momentum integration in the second term 
of the sum on the r.h.s. of Eq. (B3), obtaining the expression 



tanh(i/3eK) 



tanh(i/3e-K) 



(B5) 



Two-dimensional integration in (B5) reads as 
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dK= KdK / (B6) 
Jo ^0 

Using the previously defined dimcnsionless parameters [cf. notations ( |A^ )] along with new variables 



and the dimensionlcss function C [cf. notations (A?)] 

C(Q„0) = j.oC(q,) (B8) 

one gets 

= /%fcdfctanh/3o(fc'-l) r^ ^., , ^ a (B9) 



— / c?z tanh Bnz — ^ ^ 

4^7-1 4(z + l)g2 



C(q,) = — / dzt&nhPoz ^ la (BIO) 



We denote la according to Eq. ( |A9| ) with 

- > (BID 

4(z + l)g2 - 

Note that 

After some simple manipulations we get 

= i / " , , ''"^f f = sign(z + 29f ) e(.^ > 4<z,^(l - g?)) (B13) 

The square root is chosen to be positive on the real axis of integration z. Recalling that g2 — cos2 6 and, consequently, 

4g«'(l-9?)-sin2 20, (B14) 

we see that the sign- and event functions play their role on the patch z G [—1,1] only. Due to the symmetry properties 
of function C{qs = cos 6') it suffices to consider the angle range < 6* < 7r/2. After some analysis one can obtain 

^/ X 1 /"'^ , tanh/3o2 „/7r ^ ^ 7r\ /"^ , tanh/^o^ /-r^,,-N 
C{qs)^^ dz ^" +9/ <g< / dz ' = B15 

2^1 - sin^ 261 ^4 2/7,i„2e ^^2 - sin^ 29 

The function C{qs) is continuous. It can be easily calculated at zero temperature, giving the known result (See. e.g., 
Rcf. [|i): 



— In ^ 

2 q^s 



C(g,) = -ln— , 0<g, <1 , (B16) 



as well as for zero total incoming momentum {qs — 0) at low (To ^ 1) temperature, giving with exponential precision 

C(0)=ln(l^A) (B17) 

where 7 = 0.577.. is Euler's constant. For the latter calculation cf., e.g., section 33.3 of Ref. [p. In the above formulas 
and in what follows we track A up to terms 0(1). 
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2. Fourier components 

From Eq. ( B15| ) it is easy to calculate Fourier components of function C'{qs) given by 



C„ = - [\e cos 2n0 c{e) 



For the zeroth component after some manipulations we have 



1 r 

Cq = — dztanh/?oz 

Jo JO 

The first term can be written as 



\J — sin^ V "'0 



tanh (5^z 
\/ — sin^ V 



1 /"^ 

Co^^ = - / d^tanh/3oz 

where K(z) is the complete elliptic integral of the first kind, which has the following Taylor expansion: 



To proceed further we notice that 



{2' 2' 2}'{^' 2}'' 



F(z) 



(B18) 



(B19) 



(B20) 



(B21) 



(B22) 



where the Taylor expansion of the hypergeometric function can be readily read from the term by term integration of 
the expansion (B21), i.e., 



^ ' 2 V 12 320 ^ V 



Integrating by parts the r.h.s. of Eq. (B20) we obtain 



Ci^) = i(F(l).tanh/3o- /' ^"f^ F(z)) 
TT V Jo cosh fi^z ' 



(B23) 



(B24) 



Using the same (exponentially accurate at low temperature) approximations which lead us to Eq. (A26), we end up 

dz /„ Tl 3 grs 



TT z ./n cosh 



(1) 



(B25) 



In order to calculate F(l) it is convenient to come back to the integral representation of the elliptic integral, and then 
to evaluate it via changing the order of integration. 



, , r r dx 1 
F(l)= / dz \ -== 
Jo Jo \J \ — \J z'^ — 

Here G — 0.91569... is Catalan's constant. Finally 

2 

TT 



dx 



dz 



In 2. 
2" 



Vl — Ja; V-z^ — 



3C(3), 



= 2 / de'lncot6l = 2G 



For the second contribution Cq one integration by parts gives 

1 /"^.V. 2A2 /■^'"Vorfzln(z+ v/z2^^^h?V')^ 



- / ' dv/ln- 

TT ./n V 1 



COS-0 



cosh f5oz 



(B26) 



(B27) 



(B28) 



It is easy to see that the second term on the r.h.s of Eq. ( B28| ) is exponentially small, while the first one can be readily 
evaluated resulting in 
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Ci^^ = ln2A- -(! 



(B29) 



Combining our results together, we have 

Co = ln2A 



In 2 



-Tn 



3C(3) 
64 



Calculation of the first Fourier component Ci is very simple, since the first term on the r.h.s. of Eq. (BK 
after the angular integration, while the second one can be evaluated exactly, resulting in 

Ci = lncosh/?o = + + 0(e-i/^°) 



(B30) 
gives zero 

(B31) 
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FIG. 1. Diagrammatic equation for the four-point vertex at one-loop level. The one-loop graphs are called ZS, ZS', and 
BCS in the order they appear on the r.h.s. of this equation. 
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